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^ ! Abstract 

O 



Following the subtraction procedure for manifolds with boundaries, we cal- 
culate by variational methods, the Schwarzschild and Flat space energy dif- 
ference. The one loop approximation for TT tensors is considered here. An 
analogy between the computed energy difference in momentum space and the 

X" 

;h ' Casimir effect is illustrated. We find a singular behaviour in the UV-limit, 

due to the presence of the horizon when r = 2m. When r > 2m this singu- 
lar behaviour disappears, which is in agreement with various other models 
previously presented. 

I. INTRODUCTION 

An interesting problem appearing in Einstein gravity is the computation of quantum 
corrections to a classical energy. A possible approach is the analysis of the thermodynamical 
quantities that characterize the system under consideration. This analysis can be carried 



out by computing the system free energy at a given volume and temperature, by means 
of a partition function and the Euchdean action. Following the background method, we 
fix a metric and look at quantum fiuctuations with respect to such a background with the 
appropriate boundary conditions, then we functionally integrate such metric fiuctuations 
which are strictly periodic in Euclidean time t. In particular, the only feasible way to treat 
functional integration is by saddle-point methods. This is adequate for the treatment of the 
small perturbation concerning Minkowski space and for a semiclassical analysis of vacuum 
stability. However, a different point of view based on the Hamiltonian approach, could be 
considered. In this framework, quantum corrections to a classical energy can be computed 
by means of expectation values of the total Hamiltonian with respect to some states. It is 
clear that the problem is too large to be completely solved. To this end we might take into 
consideration the simplest non trivial saddle-point we can extract from vacuum Einstein 
equations, the Schwarzschild solution 

2MG\ ,, / 2MG\-^ 



ds' 



1 \dt^+\l j dr^ + Pdn^, (1) 



where dQ"^ = d9'^ + sin^ 9d(/)'^ is the line element of the unit sphere, G is Newton's constant 
and M is a parameter representing the mass of the wormhole. This metric is asymptotically 
fiat. The apparent singularity located at r = 2MG can be removed by a suitable definition 
of the coordinates, e.g., the Kruskal-Szekeres coordinates, which is written as 



2MG ^)'^^'P[2Mg) ~ ^y (2) 

^^P V2MG j ^ y 

In terms of these coordinates we have 

ds' = ^li^ exp (-^) dxdy + r'dn\ (3) 

The only true singularities are at curves xy = —1, where r = 0. The region {x > 0,y > 0} 
is the ^^ outside region" , the only region from which distant observers can obtain any infor- 
mation. The line y = 0, where r = 2MG, is the ^^ future horizon"; the line x = where also 
r = 2MG, is the ^^past horizon" . We will consider a slice S of the Schwarzschild manifold 



representing a constant time section of A4. This surface S is an Einstein- Rosen bridge with 
wormhole topology 5^ x R^ which defines a bifurcation surface, dividing S in two parts de- 
noted by S+ and S_. Our purpose is to consider perturbations at S with t constant, which 
naturally define quantum fluctuations of the Einstein-Rosen bridge. In particular we will 
focus our attention on the E+ sector of the manifold, corresponding to the ^^outside region" 
of the Kruskal manifold. The explicit expression of the Hamiltonian can be calculated by 
considering the line element 

ds'^ = -N'^ [dx^Y + gij [N'dx^ + rfx^) (N^dx^ + dx^) , (4) 



where A^ is called the lapse function and Ni is the shift function. When A^ 



2MG 



Ni = and gijdx^dx^ = ( 1 — ^^^) dr"^ + r'^dVt^^ we recover the Schwarzschild solution. 
On the slice S, deviations from the Schwarzschild metric spatial section will be considered 

9ij = 9ij ~\~ hij, (p) 

with A^j = and N = N (r). Then the line element (^) becomes 

ds"^ = -N^ (r) (^dx'^y + gijdx'dx^ (6) 

and the total Hamiltonian is 

Ht = H^ + Haj, = ^ d^'xiNn+N.W) + Rgj,, (7) 



where 



l-i=G,jki^'^Tx^^ ( ^ ) - ( ^ ) i?^^^ (Super Hamiltonian), (8) 



.V9) V ^: 



rC = -27r[j (Super Momentum), (9) 

while if as represents the energy stored in the boundaries. In this respect, we will follow the 
AVAi approach ||^, even though the quasilocal energy context gives a more general treat- 
ment with the possibility of looking at the gravitational thermodynamics MM- Moreover, 



since the space under investigation is asymptotically flat in spacelike directions, the quasilo- 
cal energy agrees with the results of the AVM. approach in the limit that the boundary 
tends to spatial infinity. In any case, to correctly compute Hqy, = Hj^^j^ we have to fix 
a reference frame to normalize the energy value on the boundary to zero. This opens the 
problem of the subtraction procedure investigated in Refs. |§,§]. In this paper we would 
like to apply such a procedure extended to the volume term, at least at one loop. Since 
the reference space for the Schwarzschild metric is fiat space, the contribution to the energy 
term is 

Havm = lim f Jdg"^ [9iKj - Qtj.k] dS^ = M (10) 

where gij is the metric induced on a spacelike hypersurface (9S which has a boundary at 
infinity like S^. Following the result of Ref. 0, we see that Hj^j^m is completely equivalent 
to 



' ' ■'" ■'" (11) 



2t^ 2 



K-'K, 



8nG JdT. 

where the subtraction structure is evident. The one-loop contribution to the zero point 
energy for gravitons embedded in fiat space is 

2 ■ - / ^VF. (12) 

2i (27r)' ^ ^ 

It is clear that this term is UV divergent. We will show that the same kind of divergence is 
present when the curved background is considered. In the spirit of the subtraction procedure 
we will compute the difference between zero point energies. Their difference at one loop 
represents a Casimir-like computation. The paper is structured as follows, in section |I| we 
define the gaussian wave functional for gravity and we analyze the orthogonal decomposition 
of the metric deformations, in section |IT1| we give some of the basic rules to perform the 
functional integration and we define the Hamiltonian approximated up to second order, in 
section |V|, we analyze the spin-two operator acting on transverse traceless tensors, only for 
positive values of E"^ . We summarize and conclude in section 0. 



II. ENERGY DENSITY CALCULATION IN SCHRODINGER 

REPRESENTATION 

As already mentioned, we would like to discuss the possibility of generalizing the bound- 
ary subtraction procedure. To this end, by looking at the Hamiltonian structure, we see 
that there are two classical constraints 

n =0 

(13) 

which are satisfied both by the Schwarzschild and Flat metric and two quantum constraints 

7^^ =0 

(14) 

W^ = 

7i\E' = is known as the Wheeler-DeWitt equation (WDW). Nevertheless, we are interested 
in assigning a meaning to 



^ 



zjSchw. zjFlat 



^) , {^\H^vm\^)^ (^5) 



where \& is a wave functional whose structure will be determined later and Hj^^'^'^' (hP""^) 
is the total Hamiltonian referred to the different spacetimes for the volume term. This 
has to be meant in this way: it is true that the WDW equation refers to the space of 
metrics, but the space of metrics posses different sectors [|l[ and we are considering the 
sector of asymptotically flat metrics, in which the zero point energy is deflned with respect 
to Minkowski space. For the de Sitter sector, we have to substract the energy of de Sitter 
background and so on. Note that if the expectation value is calculated on the wave functional 
solution of the WDW equation, we obtain only the boundary contribution. However, in this 
context boundaries are at inflnity in spacelike directions, that it is equivalent to considering 
the unphysical situation of computing energy excitation in the asymptotic region. Then 
to give meaning to (P^, we adopt the semiclassical strategy of the WKB expansion. By 
observing that the kinetic part of the Super Hamiltonian is quadratic in the momenta, we 
expand the three-scalar curvature J (Pxy/gR^^^ up to o{h?) and we get 

5 



/ 



d X 



-^hAh + h'^Ah,, - i/i*^ VzV,/i5 + ^hViV^h'^ - ^h^'Riah"; + \hR,,h'^ 



(16) 



where h is the trace of hij. On the other hand, following the usual WKB expansion, we 
will consider ^ ~ Cexp(zS'). In this context, the approximated wave functional will be 
substituted by a trial wave functional according to the variational approach we would like 
to implement as regards this problem. 

III. THE GAUSSIAN WAVE FUNCTIONAL FOR TT TENSORS 

To actually make such calculations, we need an orthogonal decomposition for both tTj^ and 
/ijj to disentangle gauge modes from physical deformations. We define the inner product 

{h, k) ■= f ^G'^^'h^j (x) hi (x) d^x, (17) 

Jm 

by means of the inverse WDW metric Gijki, to have a metric on the space of deformations, 
i.e. a quadratic form on the tangent space at h, with 

Qijkl ^ ^gikgjl _^ gilgjk _ 2g''Jg^^). (18) 

The inverse metric is defined on co-tangent space and it assumes the form 

(P, g) := / VgGijkip"^ (x) q^^ (x) d^x, (19) 

SO that 

G'^'^'^GnmH = \ (5151 + 5\5i) . (20) 

Note that in this scheme the "inverse metric" is actually the WDW metric defined on 
phase space. Now, we have the desired decomposition on the tangent space of 3-metric 
deformations P,|5|: 

hi, = \h9.j + {Li\j + hi^ (21) 

where the operator L maps ^j into symmetric tracefree tensors 



(LO,, = V.e, + V,6 - ^9^, (V ■ • 



(22) 



Then the inner product between three-geometries becomes 



J M 



M 



V~9 



';h' + {LO'^ (LO,, + h^^^h^ 



(23) 



With the orthogonal decomposition in hand we can define a " Vacuum Trial State" 



1 



^[hi,{lf)]=Afexp{- — 



4/2 



hK~^h 



± 



x,y 



m)K-'m) 



x,y 



hR-'h 



Trace 



x,y 



, (24) 



which will be used as a probe for the gravitational ground state. This particular expression 
is useful because the functional can be represented as a product of three functionals defined 
on the decomposed tensor field 



^ [/i,, (If)] = Af^ \hj-^ (If )] ^ \m) 



^ 



o9ijh ( x^; 



(25) 



hjj is the tracefree-transverse part of the 3D quantum field, {L^)^- is the longitudinal part 
and finally h is the trace part of the same field. {■,-)xy denotes space integration and 
K~^ is the inverse propagator containing variational parameters. The main reason for a 
similar ^^Ansat£^ comes from the observation that the quadratic part in the momenta of the 



Hamiltonian decouples in the same way of eq. (|23|) . Note that the decomposition related to 
the momenta is independent of the choice of the functional. To calculate the energy density, 
we need to know the action of some basic operators on \E' [hij]. The action of the operator 
hij on |\E') = \I' [hij] is reahzed by 



/iij(x)l^) = hij {-^) "^ [hij] . 
The action of the operator ttjj on |\E'), in general, is 

'■' '"' i*> = -'JhTm" ■''"' 

The inner product is defined by the functional integration: 



(26) 



(27) 



(^1 I ^2) = / [Vhi,] ^l {h,j} ^2 {hkl} 



(28) 



and the energy eigenstates satisfy the stationary Schrodinger equation: 

6 



d'xH { -i- 



hij (It) \ ^ {hij} = E<I/ {hij} 



(29) 



6hij ("af ) 

where Ti. I ~i §f^^.i^\ ■, h-ij (x) \ is the Hamiltonian density. Note that the previous equation 
in the general context of Einstein gravity is devoid of meaning, because of the constraints. 
However in the semiclassical context, we can give a meaning to eq. (pOf) , where a semiclassical 
time is introduced in the same manner of Refs. |T^JT^. There, a Schrodinger equation of 
the form 



dt 



H\2^^ 



(30) 



is recovered by the WDW equation approximated to second order for a perturbed minisu- 
perspace Friedmann model without boundary terms. When asymptotically flat boundary 
terms are present we have to take account of such contributions in the WKB expansion 
such as in Ref . [^ . However in this paper only gravitational transverse-traceless modes are 
considered on the fixed curved background and \E'"'" is substituted by a trial wave functional. 
To further proceed, instead of solving (P^jj), which is of course impossible, we can formulate 
the same problem by means of a variational principle. We demand that 



(31) 



(^1^) s[Vgi;\ \^{gt^}? 

be stationary against arbitrary variations of ^ {hij}. The form oi {^ \ H \'^) can be com 
puted as follows. We define normalized mean values 

S[Vgi^^S<fxgi^{-t)\^[gi^}\' 



-gi, (^) 



/ [^9i^ I ^ [ati] 



gii (^) gii Cf) + kL^ c^, if) 



j[Vg^]jd^xg^^{l^)gii{t)\^{g^,} 



8 



(32) 
(33) 
(34) 



It follows that, by defining hjj = gij — gij, we have 



Vh 



«j 



hi-^{lt)\^{hi-^ + -gl-^}\'=0 



(35) 



and 



Vh 



«j 



d'xhi-j (If) 4 it) I ^ {^^ + g^} 



^ijki [^ : y 



Vh 



«j 



^{^i + ^.^} 



(36) 



Nevertheless, the application of the variational principal on arbitrary wave functional does 
not improve the situation described by the eg. (PUJ) . To this purpose, we give to the trial 
wave functional the form 



m 



hl^=^--^[-^{^9-^)K-\g-^)]^^^,^ 



«j 



(37) 



We immediately conclude that 



^14(^)1^ 



(38) 



where vr^^ is the TT momentum. In Appendix O, we will show that 



V 



^K (^) <i (^) l*> = -.Kim (^, t) ■ 



^ijkl 



(39) 



Choice (^) is related to the form of the Hamiltonian approximated to quadratic order in 
the metric deformations. Indeed, up to this order we have a harmonic oscillator whose 



ground state has a Gaussian form. By means of decomposition (^Tj), we extract the TT 
sector contribution in the previous expression. Moreover, the functional representation (p5D 
eliminates every interaction between gauge and the other terms. Then for the TT sector 
(spin-two), one gets 



£x^R'^^^c^—jj^x^ 



h^'^ (Ao): hi - 2hRiM'' 



V' 



(40) 



where (A2)" := — A(5j + 2Rj. The latter term disappears because the gaussian integration 
does not mix the components. Then by collecting together eq.(^D|) and eq. (PD|), one obtains 
the one-loop-like Hamiltonian form for TT deformations 



The propagator K-^ (x, x)^^^i comes from a functional integration and it can be represented 
as 

where /i^ ("af ) are the eigenfunctions of Ag,- and Xn {p) are infinite variational parameters. 



IV. THE SPECTRUM OF THE SPIN-2 OPERATOR AND THE EVALUATION 
OF THE ENERGY DENSITY IN MOMENTUM SPACE 

The Spin-two operator is defined by 

(A2); := -Ad; + 2R] (43) 

where A is the curved Laplacian (Laplace-Beltrami operator) on a Schwarzschild background 
and Rj is the mixed Ricci tensor whose components are: 

where 2m = 2MG. This operator is similar to the Lichnerowicz operator provided that we 
substitute the Riemann tensor with the Ricci tensor. This is essentially due to the fact that 
the Riemann tensor in three-dimensions is a linear combination of the Ricci tensor. In ( ^5D 
the Ricci tensor acts as a potential on the space of TT tensors; for this reason we are led to 
study the following eigenvalue equation 

[-A5] + 2Rfj h\ = E% (45) 

where E"^ is the eigenvalue of the corresponding equation. In doing so, we follow Regge 
and Wheeler in analyzing the equation as modes of definite frequency, angular momentum 
and parity. We can specialize to the case with the quantum number corresponding to the 
projection of the angular momentum on the z-axis is zero, without altering the contribution 

10 



to the total energy because of the spherical symmetry of the problem. In this case, Regge- 



Wheeler decomposition |[TT| shows that the even-parity three-dimensional perturbation is 



h'r^^-{r,{},c^) = dtag 



H{r){l ) ,r^K{r),r^sm^^K{r) 



Yioi^A)- (46) 



Representation (^) shows a gravitational perturbation decoupling. For a generic value of 
the angular momentum. angular momentum /, one gets 

-AiH{r)-^H{r) = EfH{r) 



-AiK{r) + ^K{r) = EfKir) 



-AiK{r) + ^K{r) = EfK{r). 



The Laplacian in this particular geometry can be written as 

_/ 2m\ cP /2r-3m\ d 1(1 + 1) 
\ r J dr"^ \ r^ j dr r^ 



(47) 



(48) 



Defining reduced fields, such as: 



H{r) 



h (r) 



Kir] 



k (r) 



(49) 



and changing variables to 






x = 2m{ \l^\l- 1 + In {\ — + \ — - 1 

2mV2m W 2m V 2m 



(50) 



the system ( |4^) becomes 






dP 



h (x) + V- {x) h [x) = Efh (x) 



-£,k{x) + V+ix)kix) 






k [x) + V+ [x) k [x) 



Efk ix) 



Efk (x) 



where 



(51) 
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/(/ + !) 3m 



r^ (x) r (x) 



3- 



(52) 



This new variable represents the proper geodesic distance from the wormhole throat such 

that 

when r — > oo, x ~ r and V^ (x) — > 



when r — > tq, x ~ and V^ (x) 



/ (/ + 1) 3m 

2 — -F ~3" = const, 



(53) 



'0 '0 

where tq satisfies the condition tq > 2m. The solution of (^Tj), in both cases (fiat and curved 
one) is the spherical Bessel function of the first kind 



joipx) = \ -sm{px) 

IT 



This choice is dictated by the requirement that 



(54) 



h{x) ,k (x) -^ when x — i> (alternatively r -^ 2m) . 



(55) 



Then 



K{x,y) 



io (pa;) io (P?/) 1 



2A 



4tt 



(56) 



Substituting (|56|) in (]4l|) one gets (after normalization in spin space and after a rescaling of 
the fields in such a way as to absorb l^) 

V ^A /•- , of... Ef{p,m,iy 



^''^ 1=0 i=l -^0 



A.(p) 



(57) 



where 



77^2 r„ ^ 7^ _ „2 , U^ + 1) _ 3m 



rn 



T 



^3 ' 



(58) 



Aj (p) are variational parameters corresponding to the eigenvalues for a (graviton) spin-two 
particle in an external field and V is the volume of the system. 

By minimizing { ^7\ ) with respect to Aj (p) one obtains Aj (p) = [Ef {p, m, I)] ^ and 
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with 



T^ OO 2 „Q(3 

E(m,X) = —Y^Y. dp2^E!{p,m,l) 



2 / / + 1 3m 
P + 2 > — ■ 



(59) 



Thus, in presence of the curved background, we get 

V 1 °° 






(60) 



where 



2 _ / (/ + 1) 3m 



"T 



while when we refer to the flat space, in the spirit of the subtraction procedure, we have 
m = and c2 = ^^. Then 

y \ "^ roo 



Em^^-^Ef^ppi^^/^^) 



(61) 



Now, we are in position to compute the difference between ( |60D and (|6TD. Since we are 



interested in the UV limit, we have 



AE{m) = E{m) - E {Q) 



y \ ^ roo 

2^2 gio '"'' 



p^ + c± + Jp^ + c+ — 2Jp^ + c 



,2 I ^2 



y I /"oo 



\ 



-7h 



\ 



-7 - 



\ 



\p, 



and for p^ >> c^,c^. we obtain 



V 1 

2^2 



°^ POO 



^^2 U J sL J ^^^2U J 8 [p] 



(62) 



.P. 



4U. 



K c^ /"^ (ip 



27r2 8 7o p 



(63) 
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We will use a cut-off A to keep under control the UV divergence^ 



A 



^^r^^ln(^). (64) 

op Jo X \ cj 

Thus AE (m) for high momenta becomes 

.r.. ^ V c\ /An V /3mV 1 , frlA^\ 

AE(m) ^— n — = ^^ — n ^ — . (65) 

At this point we can compute the total energy, namely the classical contribution plus the 
quantum correction up to second order. Recalling the definition of asymptotic energy for 
an asymptotically flat background, such as the Schwarzschild one gets, 

M 7;\^] — n ^ — ] = M -\ — ^ — n ^-— (66) 

27r2 \r^ J 16 \ 3m J 2n^ \ r^ J 16 \3MG J ^ ' 

One can observe that 

AE (m) -^ oo when ?7i — > 0, for tq = 2m = 2GM (67) 

and 

A^ (m) -^ when m -^ 0, for tq ^ 2m = 2GM. (68) 



Remark We would like to explain the reasons that support the results of formula (pSf). In 
that formula we introduced a particular value of the radius, which behaves as a regulator 
with respect to the horizon approach of the potential. The meaning of this particular value 
is related to the necessity of explaining the dynamical origin of black hole entropy by the 
entanglement entropy mechanism and by the so-called ^^brick wall moder PJTO[]. Indeed, 



^It is known that at one-loop level Gravity is renormalizable only in flat space. In a dimensional 
regularization scheme its contribution to the action is, on shell, proportional to the Euler character 
of the manifold that is nonzero for the Schwarzschild instanton. Although in our approach we are 
working with sections of the original manifold to deal with these divergences one must introduce a 
regulator that indeed appears in the contribution of energy density. 

14 



the same mechanism is present when one has to regularize entropy by imposing a kind of 
cut-off, which in coordinate space means r^ > 2m. In fact, tq can be seen as 2m + h, where 
h assumes the same meaning of Ref. 0. However, to exphcitly relate this quantity we have 
to compare the Bekenstein-Hawking entropy with the result deriving from the evaluation of 
the free energy for gravitons, in this case. The only difference from the original calculation 
is the spin contribution not present for scalar fields. 

V. SUMMARY AND CONCLUSIONS 

We started from the problem of defining quantum corrections (semiclassical) to a grav- 
itational energy. By means of a variational approach with Gaussian wave functionals an 
attempt to calculate such a correction was made. Despite the constraint equations, this 
calculation is based on an extension of the subtraction procedure involving volume terms in 
the semiclassical regime. Excitations coming from boundary terms have been neglected to 
avoid the unphysical situation of having contributions deriving from infinity. In this context 
the extended subtraction procedure corresponds to the difference between zero point ener- 
gies calculated in an asymptotically flat curved background referring to a fiat background. 
This procedure eliminates the UV divergence of the free gravitons leaving the contribution 
of the curved background related to an imposed by hand UV cut-off. A strong analogy with 
the Casimir vacuum energy calculation is revealed, opening the possibility of understanding 
several configurations and their relationship with the vacuum stability. Indeed, this appa- 
ratus can be applied also to the Schwarzschild-deSitter background which asymptotically 
approaches the deSitter space and so on. Although this evaluation has to be completed with 
a careful study of the spin-two operator, we can conclude that the variational approach for 
the computation of quantum corrections (semiclassical) to a classical energy can be thought 
of as a tool for testing zero point energy (Casimir-like energy) in a complicated theory such 
as Einstein gravity 
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APPENDIX A: CONVENTIONS 

1. Riemann tensor, Ricci tensor and the Scalar Curvature in 3Z) 

4- = rLj - r;.,„ + rj^r^, - Vlj% Riemann tensor 

Because of the vanishing of the Weyl tensor in 3D, that is C'^^ = 0, the Riemann 
tensor is completely determined by the Ricci tensor 

Rim = Rum Ricci tensor 
R = g''-' Rij Scalar curvature 

APPENDIX B: THE KINETIC TERM 

The Schrodinger picture representation of the kinetic term is 

V 5hij (x) 5hki (x) J 
We have to apply this quantity to the gaussian wave functional |\E'). This means that 

tt'^ (x) TT^' (x) I^) ^ ^ ^ 



6hij (x) 6hki {x) 



2 
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-- fd'y'dV (^y 49¥)\[^)K~'^''^^'''^ {x,y') h,'i' {y') 



.K-'^^^^^'"''^{x,y")h,„i"{y")^[h]. (B2) 

By functional integrating 

(^ \hkn' iy') hu'n" {y")\ ^) = i^(fe'/')(fe"P) (?/', /) (^1^) • (B3) 



Thus 



\1> tt'^' ix) tt'"'- (x) ^ 



becomes 



^K-'^mm^x,x)(^) 



\ f d'y'dY ( V^) ' ^[^)^[^)K-'^''^^''^'^ {x, y') K~^^^)^'^") (x, y") 



Kik'i'){k"i"){y\y"){^\^) 

1 r.-li 



-X-i('=')(^^)(a:,x)(v^) (^1^)- (B4) 

Then the expectation value of the kinetic term, with the Planck length reinserted, is 

(T) = ^ I d\^ (G^.mK^'^"'^^'^^ (x, x)) , (B5) 
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